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All  crocodilians  are  under  various  threats  due  to  over  exploitation  and  these  species  have  been  listed  in 
Appendix  I  or  II  of  CITES.  Lack  of  molecular  techniques  for  the  forensic  identification  of  confiscated 
samples  makes  it  difficult  to  enforce  the  law.  Therefore,  we  herein  present  a  molecular  method  devel¬ 
oped  on  the  basis  on  16S  rRNA  gene  of  mitochondrial  DNA  for  identification  of  crocodile  species.  We  have 
developed  a  set  of  16S  rRNA  primers  for  PCR  based  identification  of  crocodilian  species.  These  novel 
primers  amplify  partial  16S  rRNA  sequences  of  six  crocodile  species  which  can  be  later  combined  to 
obtain  a  larger  region  (1290  bp)  of  16S  rRNA  gene.  This  16S  rRNA  gene  could  be  used  as  an  effective  tool 
for  forensic  authentication  of  crocodiles.  The  described  primers  hold  great  promise  in  forensic  identifi¬ 
cation  of  crocodile  species,  which  can  aid  in  the  effective  enforcement  of  law  and  conservation  of  these 
species. 

©  2012  Elsevier  Ltd  and  Faculty  of  Forensic  and  Legal  Medicine.  All  rights  reserved. 


1.  Introduction 

Crocodiles  are  the  largest  species  among  the  members  of  class 
Reptilia  in  vertebrates  and  are  one  of  the  last  living  links  to  the 
era  when  dinosaurs  lived.  Until  the  early  1950s,  these  species 
were  numerous  in  the  tropics  of  Africa,  the  Americas,  Asia  and 
Australia.  In  the  recent  years  many  of  crocodilian  species  are 
considered  endangered  due  to  their  habitat  destruction  and  illegal 
poaching  for  their  lucrative  products.34,35  According  to  Conven¬ 
tion  on  International  Trade  in  Endangered  Species  (CITES)  of  Wild 
Fauna  and  Flora,  crocodilian  species  were  grouped  in  Appendix  I 
and  II.  A  significant  illegal  trade  of  crocodilians  has  engrossed 
a  lot  of  attention  and  therefore  the  species  deserve  priority  for 
conservation.  In  order  to  conserve  these  antiquities  many  coun¬ 
tries  have  developed  legal  policies  for  prohibition  of  poaching  and 
hunting.28  However,  illegal  trade  still  exists  in  many  parts  of 
world  which  appears  to  be  a  challenging  task  to  proscribe. 
Moreover,  the  biological  materials  are  often  degraded  quality, 
which  do  not  lead  to  conclusive  identification  through  conven¬ 
tional  methods.6,39  In  this  perspective,  development  of  novel 
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methods  for  the  identification  of  these  species  is  inevitable 
towards  their  conservation. 

Mitochondrial  (mt)  DNA  is  valuable  tool  for  species  identifica¬ 
tion  from  highly  degraded  samples.1,11  Many  of  universal  mt  DNA 
markers  such  as  cytochrome  i>,18,38  12S  rRNA,17  16S  rRNA,14  and 
noncoding  D-Loop  region28,40  have  shown  their  potential  targets 
for  animal  species  identifications.  Within  crocodiles  a  single 
mtDNA  marker  i.e.  cytochrome  b  (cyt  b )  has  been  studied  as 
a  potential  source  for  species  identification.21  ~23, 25,27  However,  the 
efficiency  of  other  regions  for  forensic  identification  of  crocodiles 
has  not  been  evaluated  and  thus  these  markers  also  need  to  be 
studied  in  forensic  analyses. 

The  16S  rRNA  gene,  which  encodes  the  mitochondrial  large 
ribosomal  subunit  (mt  LSU)  in  animals,  has  been  employed 
extensively  to  explore  phylogenetic  relationships  in  anthropods  at 
most  phylogenetic  levels,8  familial  level2  and  the  genus  level  and 
below.5,7,15  The  wide  range  in  utility  of  16S  at  various  taxonomic 
levels  suggests  that  it  could  be  useful  marker  for  forensic  identifi¬ 
cations  as  well. 

In  this  paper,  we  develop  novel  primer  sets  to  amplify  partial 
16S  ribosomal  RNA  region  for  discriminating  the  crocodile  species. 
The  primers  were  optimized  and  checked  for  their  efficiency  in  six 
species  i.e.  Crocodylus  palustris,  Crocodylus  porosus,  Gavialis 
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gangeticus,  Crocodylus  niloticus,  Crocodylus  siamensis  and  Croc- 
odylus  crocodilus  which  represent  the  three  existing  families  of 
Crocodylia. 

2.  Materials  and  methods 

2.2.  Sample  collection  and  DNA  extraction 

Authenticated  biological  samples  were  obtained  from  Madras 
Crocodile  Bank  Trust  (MCBT),  Centre  for  Herpetology,  Mamalla- 
puram,  Tamilnadu,  India,  National  Chambal  Sanctuary  Project,  Agra, 
Uttar  Pradesh,  India,  Nehru  Zoological  Park,  Hyderabad,  Andhra 
Pradesh,  India,  under  the  consent  of  Ministry  of  Environment  and 
Forests,  Government  of  India,  New  Delhi.  Whole  blood  samples  of  C. 
palustris,  C.  porosus,  Crocodylus  siamensis,  C.  niloticus,  Caiman  croc¬ 
odilus,  and  G.  gangeticus  were  provided  by  Madras  Crocodile  Bank 
Trust  (MCBT),  Tamilnadu,  India  and  15  dead  gharial  samples  (fresh 
tissue,  putrified  tissue  and  bone)  were  provided  by  National 
Chambal  sanctuary,  Utter  Pradesh,  India.  Sample  numbers  and 
breed  of  each  species  are  listed  in  Table  1.  Total  Genomic  DNA 
extraction  from  blood  samples  were  carried  out  by  standard 
phenol-chloroform  procedure33  and  cleaned  using  Microcon  100 
centrifugal  filter  column  (Millipore).  DNA  extraction  from  tissue 
samples  was  performed  using  Qia  tissue  DNA  extraction  kit  (Qia- 
gen,  Valencia,  A)  as  per  the  manufacturer’s  guidelines,  The  DNA 
extracts  were  adjusted  to  50  pL  in  volume  with  TE~4  buffer  and  kept 
at  -80  °C  until  use.  The  quality  and  quantity  of  all  DNA  samples 
isolated  from  blood  (100—200  ng)  and  tissue  (50  ng)  was  checked 
in  1%  agarose  gel. 

2.2.  Primer  design  and  PCR  amplification 

In  order  to  construct  these  primers,  the  complete  mt  DNA 
sequence  of  available  crocodiles  species  were  retrieved  from  NCBI 
(http://www.ncbi.nlm.nih.gov)  database  (Table  2)  and  aligned 
using  Clustal  W  program36  and  MEGA  5.37  Two  of  the  highly 
conserved  regions  in  16S  rRNA  gene,  in  all  crocodiles  species  were 
chosen  and  the  primers,  Set  I:  16SFI:  5'-AAA  GCATTC  TGC  CTA  CAC 
CTG  AAA-3';  16SRI:  5'-TTG  TGT  TGG  CTG  CIT  TAA  GGC  CTA  -3'  and 
Set  II:  16SFII:  5'-TGCTTG  GGA  ACA  GAA TAT  AAG  TTC-3' ;  16SRII:  5'- 
CTG  CIT  TTG  CAC  AGG  GAG  ATC  AAT  TTC-3'  were  designed. 

PCR  reactions  were  carried  out  using  a  total  reaction  volume  of 
20  pi  containing  2.0  pL  of  dNTPs  (2.5  mM  each)  (MBI  Fermentas, 
Glen  Burnie,  MD),  2.0  pL  of  10X  buffer  (Containing  200  mM  Tris— 
HC1,  pH  8.4,  and  500  mM  KCl)  (Invitrogen  Life  Technologies),  0.6  pL 
of  Exprime  Taq  DNA  polymerase  (5  U/pL)  (Genet  Bio),  and  1.0  pL  of 
genomic  DNA  (100  pg).  Sterile  water  was  added  to  make  up  the 
volume  to  20  pi.  All  the  PCR  reactions  were  carried  out  on  Gen- 
Amp  9700  thermal  cycler  (Applied  Biosystems,  Foster  City,  CA) 
under  the  following  cycle  conditions:  94  °C  for  5  min  of  initial 
denaturation  followed  by  30  cycles  of:  denaturation  at  95  °C  for 


Table  1 

List  of  animals  investigated  in  this  study  ( N  —  No.  of  individual  samples). 


Species 

N  (total) 

Breed 

ID  number 

Gavialis  gangeticus 

10  (Blood) 
15  (Tissue) 

Indian 

G1B-G10B,  G1T-G15T 

Crocodylus  porosus 

09(Blood) 

Indian 

S1B-S9B 

Crocodylus  palustris 

11  (Blood) 
04  (Tissue) 

Indian 

M1B-M11B 

M1T-M4T 

Crocodylus  siamensis 

06  (Blood) 

Indonesia 

CSID1-CSID6 

Crocodylus  niloticus 

06  (Blood) 

African 

CNIA1-CNIA6 

Caiman  crocodilus 

06  (Blood) 

South  American 

CCSA1-CCSA6 

Table  2 

List  of  crocodile  species  and  their  accession  numbers  as  retrieved  from  NCBI 
database. 


Species  name 

Accession  number 

Tomistoma  schlegelii 

AJ81045 

Osteolaemus  tetraspis 

AM493868 

Alligator  mississippiensis 

Y13113 

Alligator  sinensis 

AF5 11507 

Melanosuchus  niger 

EU161678 

Paleosuchus  palpebrosus 

AM493870 

Paleosuchus  trigonatus 

AM493869 

30  s;  annealing  at  63  °C  (Set  I  primers),  66  °C  (Set  II  primers)  for 
30  s,  extension  at  72  °C  for  30  s  and  amplification  ended  with 
a  7  min  final  extension  step  followed  by  a  4  °C  hold.  The  amplicons 
were  checked  in  2%  agarose  gel  along  with  a  negative  control. 

2.3.  DNA  sequencing  and  data  analyses 

The  PCR  products  were  sequenced  using  BIG-DYE  terminator 
cycle  sequencing  kit  (Applied  Biosystems,  Foster  City,  CA)  using 
the  both  forward  and  reverse  primers.  Sequencing  was  carried  out 
on  Applied  Biosystems  (Foster  City,  CA,  USA)  3100  Avant  Genetic 
Analyzer  according  to  manufacturer’s  protocols.  The  sequences  of 
two  amplicons  were  edited  using  Bioedit  software12  and  the 
sequences  were  then  used  for  phylogenetic  analysis.  The  molec¬ 
ular  distances  were  calculated  using  MEGA  5  software  with  all 
available  nucleotide  substitution  models.  Two  different  types  of 
phylogenetic  analyses,  maximum  likelihood  (ML)  and  Bayesian 
analysis,  were  performed  using  PhyML9  and  MrBayes  Version  3.1, 13 
respectively. 

3.  Results  and  discussion 

Forensic  laboratories  often  receive  degraded  samples  which  fail 
to  yield  longer  sequences,11  but  long  read  sequences  can  be  highly 
beneficial  for  the  flawless  identification  of  species.  During  molec¬ 
ular  study  of  degraded  samples  based  on  mitochondrial  DNA 
(mtDNA)  the  cytochrome  b  gene  has  been  found  to  be  a  powerful 
marker  for  identification  of  species21,31,41  and  also  used  in  studies  of 
molecular  evolution.18,29  Conversely  the  level  of  polymorphism  was 
high  in  cytochrome  b  gene  compare  to  16S  rRNA  gene,  but  species 
specific  motifs  were  higher  in  16s  rRNA  gene.10  Therefore,  we  have 
designed  two  sets  of  novel  primers  within  16S  rRNA  gene  region  of 
mitochondria.  The  two  primer  pairs  (set  I  and  set  II),  proposed 
herein  amplified  the  partial  16S  rRNA  gene  sequences  (600bp  and 


Fig.  1.  Two  percent  agarose  gel  showing  PCR  products  of  six  crocodile  species,  Lanes 
1-6  represent  600bp  amplicons  obtained  using  primers  16S  FI  and  16S  Rl.  Lanes  8-13 
represent  690  bp  amplicons  obtained  using  primers  16S  F2  and  16S  R2.  Lane  M  -  1  kb 
molecular  size  standards.  Lanes  1  and  8  —  G.  gangeticus,  Lanes  2  and  9  —  C.  porosus, 
Lanes  3  and  10  -  C.  palustris,  Lanes  4  and  11  -  C.  niloticus,  Lanes  5  and  12  — 
C.  crocodilus,  Lanes  6  and  13  —  C.  siamensis,  Lanes  7and  14  —  respective  negative 
controls. 
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Table  3 

List  of  the  intraspecies  polymorphic  sites  in  six  species. 


Species  (n  in  total) 

Position 

Nucleotide 

n 

Gavialis  gangeticus  (25) 

10a 

T 

22 

A 

3 

31 

G 

20 

A 

5 

443 

T 

24 

A 

1 

983 

G 

21 

T 

4 

1035 

G 

23 

C 

2 

Crocodylus  porosus  (09) 

43 

T 

7 

C 

2 

222 

G 

8 

A 

1 

346 

T 

8 

A 

1 

403 

A 

7 

G 

2 

637 

A 

6 

G 

3 

852 

C 

8 

T 

1 

Crocodylus  palustris  (15) 

390 

G 

14 

A 

1 

660 

G 

12 

A 

3 

888 

C 

11 

A 

4 

993 

G 

12 

T 

3 

Caiman  Crocodilus  (06) 

41 

A 

4 

C 

2 

124 

T 

5 

C 

1 

223 

T 

4 

A 

2 

908 

G 

5 

C 

1 

Crocodylus  siamensis  (06) 

31 

C 

4 

A 

2 

105 

A 

4 

G 

2 

211 

T 

5 

C 

1 

702 

A 

3 

G 

3 

Crocodylus  niloticus  (06) 

9 

G 

5 

A 

1 

359 

G 

4 

T 

2 

1029 

G 

3 

A 

3 

3  Nucleotide  numbering  started  from  the  first  nucleotide  next  to  the  3'-terminal 
of  forward  primer  in  all  species. 


690bp  respectively)  from  all  six  crocodile  species  (Fig.  1).  These 
sequences  were  later  combined  to  obtain  1290bp  sequence  of  16S 
rRNA  gene  in  six  crocodile  species  i.e.  C.  gangeticus,  C.  porosus, 
C.  palustris,  C.  siamensis,  C.  niloticus  and  C.  crocodilus.  The  sequence 
alignment  results  show  a  certain  number  of  polymorphic  sites  in 
each  species  (Table  3),  intraspecific  polymorphisms  are  an  impor¬ 
tant  deliberation  for  species  identification.40  To  validate  the 
sequence  reliability  within  each  species,  nucleotide  sequences  of 
partial  region  16S  rRNA  within  1290bp  were  analyzed  in  multiple 
individuals  of  each  of  the  six  species,  for  simplicity,  nucleotide 
position  (np)  was  numbered  by  calling  a  site  next  to  3'-terminal  of 
the  forward  primer  as  number  1  in  all  species.  A  certain  number  of 
polymorphic  sites  were  found  in  every  species,  and  three  of  the 
species  i.e.  Crocodylus  palustris,  Caiman  Crocodylus  and  Crocodylus 
siamensis  showed  less  than  four  polymorphic  sites,  similarly  G. 
gangeticus  showed  five  polymorphic  sites,  C.  porosus  showed  six 
polymorphic  sites  and  C.  niloticus  has  shown  three  polymorphic 
sites.  Most  of  the  polymorphic  sites  were  observed  in  low 
frequency.  Additionally  we  studied  nucleotide  diversity  by  using 
partial  16S  rRNA  gene  sequences  among  six  crocodile  species  are 
summarized  (Table  4).  The  percentage  range  of  sequence  diversity 
was  6.0%— 40.7%,  and  overall  mean  diversity  24.4%.  We  compared 
these  values  with  another  forensic  marker  i.e.  Cyt  b  nucleotide 
diversity  which  was  2.4%— 22.3%  and  overall  mean  diversity  20.1%. 
These  outcome  apparently  show  evaluate  to  Cyt  b  gene  which  is 
commonly  used  for  species  identification, 30,31,42  16S  rRNA  is  also 
an  important  marker  for  discriminating  crocodile  species.  The 
sequence  data  was  further  used  for  construction  of  phylogeny. 

Phylogenetic  analysis  may  help  in  resolving  the  problems  per¬ 
taining  to  sequence  consistency  and  thus  may  provide  an  additional 
accuracy  check  on  the  obtained  sequences.  Crocodilian  phylogeny 
has  been  already  examined  in  many  studies,  which  have  helped  to 
establish  higher  level  phylogenetics  relationships  between  croco¬ 
dile  species.3,4,16,19,20,24,26,32  Short  sequences  are  not  often  useful  for 
studying  the  deep  evolutionary  divergence  with  in  species14; 
therefore  we  utilized  our  newly  obtained  1290bp  16S  rRNA 
sequences  for  the  analyses  and  were  aligned  with  available 
sequences  of  other  crocodile  species  in  Genbank  using  MEGA  5.37 
Our  phylogeny  results  obtained  from  the  partial  16S  rRNA  region 
analyses  supports  the  monophyly  of  Indopacific  crocodilians  (Fig.  2) 
as  reported  by  previous  studies.23,25  The  phylogenetic  tree  obtained 
in  this  study  did  not  deviate  from  the  previously  established 
crocodilian  phylogenetic  relationships,24  which  in  turn  verify  the 
accuracy  of  the  sequences  obtained  by  combining  two  short 
sequences  of  16S  rRNA  gene.  This  also  adds  a  further  advantage  on 
the  utility  of  these  sequences  in  forensic  analyses.  Therefore  these 
novel  primers  can  be  used  to  obtain  partial  16S  rRNA  gene 
sequences  from  degraded  confiscated  samples  derived  from  croc¬ 
odile  species.  The  result  of  this  study  shows  potential  features  of 
these  16S  rRNA  primers  in  forensic  analysis  for  crocodiles.  Flence, 
we  recommend  using  these  16s  rRNA  primers  in  the  forensic 
identification  of  endangered  crocodiles  species,  which  will  help  in 
the  effective  implementation  of  wildlife  regulations  to  conserve 
these  keystone  species  from  extinction. 


Table  4 

Nucleotide  sequence  diversities  in  the  partial  region  of  16S  rRNA  (below  diagonal)  and  Cyt  b  (above  the  diagonal)3  between  six  crocodile  species  examined  in  this  study. 


G.  gangeticus 

C.  porosus 

C.  paulstris 

C.  niloticus 

C.  crocodylus 

C.  siamensis 

G.  gangeticus 

16.7% 

16.5% 

17% 

19.8% 

16.7% 

C.  porosus 

24.2% 

6% 

7% 

22.3% 

2% 

C.  paulstris 

20.6% 

6% 

7% 

20.5% 

7.5% 

C.  niloticus 

23% 

7.8% 

7% 

22.2% 

8% 

C.  crocodylus 

40.7% 

22.3% 

20.5% 

22.2% 

22.3% 

C.  siamensis 

22.6% 

12.3% 

7.8% 

8% 

22.3% 

3  Cyt  b  Nucleotide  sequence  data  was  collected  from  NCBI. 
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Fig.  2.  Bayesian  tree  with  Bayesian  probabilities  based  on  16S  rRNA  data  as  obtained  using  partial  sequences. 
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